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The classical Heisenberg model on the trillium and distorted windmill lattices exhibits a degen- 
erate ground state within large-A'' theory, where the degenerate wavevectors form a surface and 
line, in 3-dimensional space, respectively. We name such states partially ordered to represent the 
existence of long-range order along the direction normal to these degenerate manifolds. We inves- 
tigate the effects of thermal fluctuations using Monte Carlo (MC) methods, and find a first order 
transition to a magnetically ordered state for both cases. We further show that the ordering on the 
distorted windmill lattice is due to order by disorder, while the ground state of the trillium lattice 
is unique. Despite these different routes to the realization of low temperature ordered phases, the 
static structure factors obtained by large-A'' theory and MC simulations for each lattice show quanti- 
tative agreement in the cooperative paramagnetic regime at finite temperatures. This suggests that 
a remnant of the characteristic angle-dependent spin correlations of partial order remains above the 
transition temperatures for both lattices. The possible relevance of these results to /3-Mn, CelrSi, 
and MnSi is discussed. 



PACS numbers: 75.10.Hk,75.50.Ee,75.40.Cx,75.40.Mg 
I. INTRODUCTION 

Geometrically frustrated (GF) magnetic systems of- 
fer a rich avenue to the search for emergent or col- 
lective phases of matter. Frustrated magnetism arises 
when spins between nearby magnetic sites cannot form 
a unique alignment to minimize their magnetic interac- 
tions, and leads to the possibility of a macroscopically 
degenerate classical ground state. Materials with antifer- 
romagnetic nearest neighbor interactions between local 
moments lying at the sites of corner-sharing tetrahedra 
or triangles comprise a considerable fraction of the GF 
magnetic systems known to date. 

Most magnetic materials select a particular ordering 
wavevector, representative of their magnetic structure, 
upon the onset of magnetic order. In contrast, in frus- 
trated magnets, it is not uncommon to find within the 
mean field approximation a set of degenerate continu- 
ously connected "ordering wavevectors" which span the 
three-dimensional (3-D) space. For example, the disor- 
dered ground state of the classical Heisenberg model on 
the pyrochlore lattice^'^ features such a wavevector man- 
ifold. Partially ordered magnets lie between these two 
extremes, forming degenerate wavevector (2-D) surfaces 
or (1-D) lines in a 3-D space. The name of "partially 
ordered" comes from the fact that there is a long-range 
order along the direction normal to this degenerate man- 
ifold. In other words, in such systems, to very low tem- 
peratures, a system may appear disordered/ordered to 
measurements probing parallel/perpendicular directions 
to this line or surface. A growing number of frustrated 
spin models have been shown to exhibit such a partial 
order.-^i^i^ Experimental evidence for partial order has 
been observed in single-crystal neutron scattering mea- 
surements of the correlated metals, MnSi^ under pres- 



sure, and CeCus.gAuo.i^, although in these materials it 
is unlikely to arise from considerations as simple as here 
considered. 

Our interest in this problem arose from the realization 
that partial order has been found using large-A^ theory on 
both the trillium and distorted windmill lattices. The key 
question we address in this paper is whether these par- 
tially ordered states can be captured beyond the large-A^ 
theory.^ If not, to what extent can the results of large- 
A^ theory be regained due to thermal fluctuations? We 
answer this question by carrying out large scale classical 
Monte Carlo simulations of the AF Heisenberg model on 
the trillium and distorted windmill lattices. 

Large-A^ theory and MC simulations of the classical 
Heisenberg model give consistent results for both highly 
frustrated systems (ex: on the pyrochlore lattice a dis- 
ordered ground state is found), and for unfrustrated lat- 
tices, (ex: choosing the correct ordering wavevector) 
Despite this trend, the finding of a partially ordered 
ground state in large- A^ theory does not necessarily trans- 
late into a partially ordered ground state within MC sim- 
ulations. We find that the classical Heisenberg model on 
both the trillium and distorted windmill lattices shows a 
first order phase transition to a coplanar magnetically or- 
dered state featuring neighboring spins rotated by 120°. 
We further show that this model on the trillium lattice 
does not have a macroscopic ground state degeneracy. 
The chosen 120° coplanar state is a unique ground state, 
and the partial order is an artificial effect of the spherical 
approximation constraint^ used in the large- A^ theorj*^. 
On the other hand, the selection of a particular coplanar 
ordering on the distorted windmill lattice proceeds due 
to an order by disorder mechanism. 

In the cooperative paramagnetic state (for an inter- 
mediate temperature regime < T < 0cw extending 
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FIG. 1: (Color online) Magnetic sites of (a) the trillium lattice 
and (b) the distorted windmill lattice. In both cases three 
corner-shared equilateral triangles meet at each spin. The 
unit cell of the trillium lattice is only 4 sites and the next 
smallest spin loop is 5 sided. The unit cell of the distorted 
windmill lattice is 12 sites and the smallest spin loop is 4 sided. 
Atoms within the first unit cell are depicted^^ in boldface 
(black) and labeled ((b) as in Table I). 



an order of magnitude below the intercept of a Curie- 
Weiss fit, Ocw), we find that the large- description 
remains quantitatively valid for both lattices, despite the 
inability of large- iV theory to capture the true nature 
of the classical ground state. That is, calculations of 
the angle-resolved static structure factor by large- iV and 
MC techniques show at most a 5% difference. We will 
henceforth refer to this temperature window as the co- 
operative paramagnetic regime as it features reasonably 



strong but short range spin correlations, finite tempera- 
ture remnants of an avoided partial order. 

In the next section, we present the structures of the 
two lattices, and discuss the approaches used to study 
the magnetic properties of the classical Heisenberg model 
on these lattices. We show the results of MC simulations 
on the trillium and distorted windmill lattices in Sec. Ill 
and IV, respectively. In particular, we focus on the ques- 
tions addressed above: "is there a transition to a mag- 
netically ordered state at finite temperatures" and "if so, 
do we find remnants of the partial order above the transi- 
tion temperature as a result of thermal fluctuations" . We 
also address the mechanism of ordering in these lattices 
after an enumeration of their ground states. The pos- 
sible relevance of these results to real materials and the 
conclusion are discussed in Sec. V. A heuristic derivation 
of the degeneracy of the model on the distorted windmill 
lattice is presented in Appendix A. 

II. MODEL, LATTICE STRUCTURES, AND 
APPROACHES 

We study the classical 0{N) model on the trillium and 
distorted windmill lattices given by the following Hamil- 
tonian 

H = jJ2^,-S,, (1) 

where J > is the antiferromagnetic exchange coupling 
constant, the sum runs over the nearest neighbors only, 
and S = {S^, . . . , S^) is an A^-component classical spin. 

It was shown in Ref. H that the magnetic lattice 
of MnSi forms a three-dimensional network of corner- 
sharing equilateral triangles with the cubic P2i3 symme- 
try shown in Fig. 1 (a), which we have named the tril- 
lium lattice. The trillium lattice is common to many sys- 
tems including the antiferromagnetically correlated Ce 
local moments of CcIrSi. The magnetic lattice of /3-Mn, 
the distorted windmill lattice (P4i32 symmetry), bears 
a remarkable qualitative resemblance to the trillium lat- 
tice as shown in Fig. 1 (b). In particular, both lattice 
structures feature three corner-shared equilateral trian- 
gles joined at a common site. The coordinates of each 
site within a unit cell for both the trillium lattice"^ and 
the distorted windmill latticei^ have been previously pre- 
sented. For completeness, these are listed in Table H] 

To carry out a quantitative comparison between large- 
N theory and = 3 Monte Carlo results, we have car- 
ried out large scale classical Monte Carlo simulations of 
the A^=3 Heisenberg model on the trillium and distorted 
windmill lattices for lattices with Ug — nxLxLxL 
spins, where n is the number of sites (or spins) in the 
unit cell and L is the number of spins along each dimen- 
sion of a cube. On the trillium lattice (n=4), we have 
considered L = {6,9, 12, 18}. On the distorted windmill 
lattice (n = 12), we have considered L = {6, 8, 12}. The 
standard Metropolis algorithm has been used in which we 
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TABLE I: The coordinates of the trillium and distorted wind- 
mill lattice sites in the first unit cell as shown in Fig. f , are 
expressed in terms of lattice parameters u, and y respectively, 
where we set the corresponding lattice constants of the cubic 
lattices to unity. 



attempt to update a spin within a small angular range 
S around its original direction. We choose S in such a 
way that around 50% of attempted spin updates are ac- 
cepted. Starting with a random configuration, we usually 
perform 2-10^ Monte Carlo steps (MCS) for equilibration 
and 10^ MCS for measurements with one MCS consisting 
of ~ ris/T single spin updates, where T is the tempera- 
ture. 

Within large— TV theory, one writes the partition func- 
tion for the TV-component spin model— and solves for the 
spin-spin correlations'^'^^'^^ in the limit N —> oo. That is 
we rewrite the Hamiltonian as 



T 



(2) 



where Mij is the interaction matrix and the spins are 



subject to the constraint S| 
partition function is given by 



TV. The corresponding 



Z = / V(j)D\ e 



-S(0,A) 



with the action 



-M,j4>i ■ (t)j + -^5,j{4>, ■(l)i-N) 



(3) 



, (4) 



where = {(j)} , . . . , (l>^ ) is an A'^-component real vector 
field and the Lagrange multiplier for the constraint 
■ = N. To proceed, we take the oo limit 

and set a uniform = Aq. The locations i — (l,fJ-) 



of spins can be labeled by those of the cubic unit cell 
I = 1, . . . , Tic and the lattice sites /i = 1, . . . , n within the 
unit cell (ric ^ L x L x L is the total number of the unit 
cells in the lattice). The Fourier transform with respect 
to the positions of the unit cells leads to 



(5) 



with A!^" 



S^vXq- Performing Gaussian inte- 
grations over the fields in Eq. [3l one finds that Ao is 
determined by the saddle point equation. 
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(6) 



where f3e^ are the eigenvalues of the nxn interaction ma- 
trix M^" (which has been shifted such that Min(eg) = 0) 
and f3 the inverse temperature. Note that the sum over 
q is carried out for finite size periodic lattices to effec- 
tively compare with the Monte Carlo results. Plugging 
Ao back into Eqs.[3]and[51 we readily deduce the spin-spin 
correlation functions. For example, these can be found 
by calculating the second derivative of Z with respect to 
an auxiliary field which couples to the spin and can be 
added to Eq. [4] The static structure factor is found to 
be 



E 

k.k' ,p- 



(7) 



where U is the matrix that diagonalizes the interaction 
matrix M^", k and p are the sublattice indices. Note 

that, for simplicity, we will normalize spins to ^/N in the 
large- iV theoryi^i To compare with MC results (TV = 3), 
one needs to fix the energy scale. Here we rescale J 
J/3 in eP. 



III. THE TRILLIUM LATTICE 

Within large- theory, we reported^, that the classical 
antiferromagnetic (AF) Heisenberg model on the trillium 
lattice has a partially ordered ground state, with a sur- 
face of degenerate wavevectors following 



(l)-«^(f) 



where we set the lattice constant a = 1. However, 
Ref. d we relaxed the hard spin constraint of Sf = 1 



(8) 

in 
to 



where n is the number 



a soft constraint, J27=i 
of sites in the unit cell. Using energy minimization on 
spin clusters of classical spins (with A^=3), we effectively 
imposed the hard spin constraint show that the lowest 
energy state actually exhibits a coplanar magnetic order 
with wavevector (^,0,0)fi^ . We speculated that the soft 
constraint of conventional large- A?^ theory is crucial to the 
realization of partial order as T — > 0^. 
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FIG. 2: (Color online) Average energy per spin (a), heat ca- 
pacity (b) and magnetic order parameter (c) vs. temperature 
for the Heisenberg model on the trillium lattice. These results 
are independent of the lattice parameter u which determines 
the position of the sites in a cubic unit cell. In this and all 
other plots, error bars are smaller than the symbol size if not 
visible. 



To find out whether a finite temperature transition to 
a magnetically ordered state occurs, we have carried out 
large scale MC simulations. The energy as a function 
of temperature, as seen in Fig. [2](a), exhibits a smooth 
behavior before an abrupt jump, which indicates a first 
order transition to an ordered state. The heat capacity 
as a function of temperature, shown in Fig. [21(b), clearly 
shows a peak corresponding to the onset of magnetic or- 
der on the trillium lattice. Because of strong hysteresis 
effects, the location of this jump appears to scale to lower 
temperatures as the size of the system increases. We have 
not attempted to determine the precise location of this 
strongly first order transition. Our best estimate for the 
transition temperature is Tc — 0.21(1) J. In Fig. [5] (c), 
we plot the magnetic order parameter of this transition, 
which is defined as the structure factor at the ordering 



5(q = Q) 



(9) 



We see that this order parameter exhibits a sharp onset 
at the transition temperature as would be expected for a 
first order phase transition. 

It was anticipated that the MC would find a transition 
to a 120° coplanar magnetic state with the wavevector 
(^,0,0), where the spins on each triangle form at 120° 
to each other, since it was shown (by minimization) that 
this is the lowest energy state. Thus this ground state 
does not have a macroscopic degeneracy^ and the ob- 
served transition is not expected to result from an order 
by disorder mechanism. An intuitive way to understand 
the uniqueness of this ground state is to start with a 
coplanar state with spins labeled a, [3 and 7 where the 
letters a, (3 and 7 denote 120° rotated spins on every 
triangle. It is then natural to ask whether one can gen- 
erate degenerate states by rotating /3 and 7 spins around 
the axis of spin a. The high connectivity of this lat- 
tice in comparison with the kagome and hyperkagome 
lattices2S where such degeneracies naturally arise, pre- 
vents any such rotations which are not consistent with 
the crystal symmetries of the lattice. Therefore, the rota- 
tion of (3 and 7 spins around a spin axis cannot generate 
distinctly different states. This argument itself is not suf- 
ficient to prove that the ground state is unique because, 
in principle, there might be other ground states that are 
not connected to the (^,0,0) state by simple spin rota- 
tions. However, our previous minimization showed that 
the state with the wavevector (^, 0, 0) is the lowest en- 
ergy state. Combining our current MC result, we con- 
clude that the ground state of Heisenberg model on the 
trillium lattice is unique. 

In the disordered state at temperatures above this or- 
dering transition, however, the spins fluctuate, and one 
might expect to recover features found by large- A'^ theory 
with the soft constraint. Now the question is to what ex- 
tent is the partial order of large-N theory recovered due 
to temperature fluctuations? 

Magnetic susceptibility (not shown) has been calcu- 
lated within both approaches and shows good agreement 
until very close to the ordering temperature Tc « 0.21J, 
well below the Curie- Weiss temperature 0cw — 2 J, and 
in the zero temperature limit. Qualitative comparisons 
between the large- A^ theory and Monte Carlo results for 
the structure factor are shown in Fig. [3] in the coopera- 
tive paramagnetic phase at temperatures slightly above 
the onset of magnetic order. We see excellent qualitative 
agreement between these two approaches. Within the 
cooperative paramagnetic window we see that the static 
structure factor on the trillium lattice peaks around the 
surface of a sphere-like shape determined by Eq. [51 al- 
though part of this sphere is obliterated by geometric 
(rather than energetic) effects. 

Within both approaches the geometric factor (that is 
cancellations between spin contributions within the unit 
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FIG. 3: Contour plots of the intensity of the structure fac- 
tor in the hhk ((a) and (b)) and hkO ((c) and (d)) planes for 
the trillium lattice in the cooperative paramagnetic regime 
(for u = 0.138, L = 12 and T = 0.25 J) show prominant fea- 
tures near these surfaces. Classical Monte Carlo ((a) and (c)) 
agrees well with large- A'' results ((b) and (d)). The maximal 
intensity is shown in white, and the axes run from —An to An 
along k(OOl) (vertical), h(llO) (horizontal), k(OlO) (vertical), 
h(100)(horizontal). Notice the prominent features near the 
zero energy surfaces given by Eq. |8]of the large- A'^ theory. 

cell) imposes that there is very little weight within the 
first Brillouin zone. In Fig. H] a quantitative comparison 
along several high symmetry directions between these 
two approaches is presented. There is extremely good 
agreement at all wavevectors. Note that the coopera- 
tive paramagnetic regime is smoothly connected to the 
partial order within large- A'^ approaches. Therefore ex- 
cellent agreement between the angle dependent spin-spin 
correlation functions implies that remnants of partial or- 
der would be expected to appear at finite temperatures 
above the transition temperature as a result of thermal 
fluctuations. 



IV. DISTORTED WINDMILL LATTICE 

Mean field calculations'' on the distorted windmill lat- 
tice also revealed a partially ordered ground state, with a 
line of degenerate wavevectors along the {qqq) direction. 
In this paper we show that Monte Carlo simulations of 
the AF Heisenberg model on this lattice flnd an ordered 
ground state as on the trillium lattice. However, the 
physical origin of the two phase transitions differs. While 
the ordering in the trillium lattice can be understood by 
minimization, being the result of energetically different 
states, the ordering in the distorted windmill lattice can- 
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FIG. 4: (Color online) A quantitative comparison of the angle- 
dependent structure factor is shown along three high symme- 
try directions between the Monte Carlo (Heisenberg model) 
and large- A'^ results for classical spins in the cooperative para- 
magnetic phase on the trillium lattice {u = 0.138,1/ = 12). 
The solid line is the large- A*' result, while the symbols are 
from MC simulations. 



not. Rather, it must proceed via an order by disorder 
mechanism. As on the trillium lattice it is interesting 
to ask whether the partial ordering features obtained by 
large- (or mean field) theory can be found at finite T, 
where one might expect the spins to fiuctuate strongly. — 

Again, on the distorted windmill lattice, MC simula- 
tions find a jump in the energy as a function of tempera- 
ture indicative of a first order phase transition. Nonethe- 
less, this transition is qualitatively different from that 
seen on the trillium lattice, as it does not show strong 
hysteresis effects and much variation with system size. 
This is reflected in a sharp transition from one smooth 
energy vs. temperature curve to another at the transi- 
tion temperature as seen in Fig. [5] (a). The heat capac- 
ity as a function of temperature, shown in Fig. [5] (b), 
clearly shows a peak corresponding to the onset of mag- 
netic order. The variation of the order parameter shows 
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FIG. 5: (Color online) Average energy per spin (a), heat ca- 
pacity (b), and magnetic order parameter (c) vs. temperature 
for the Heisenberg model on the idealized /3-Mn lattice. 



an abrupt onset at the transition temperature. Here the 
order parameter has been defined in terms of the struc- 
ture factor at the ordering wave vector Qo — (0, 0, 0) 
as 



(m) 



2 _ 5(q - Qo) 



(10) 



As in the case of the trillium lattice, we have not at- 
tempted to determine the precise location of this first 
order transition. Our best estimate for the transition 
temperature is = 0.155(5)J. 

To understand this spin ordering we have carried out 
an energetic minimization of spin clusters with A'^ = 3. 
The minimization of spin clusters on the distorted wind- 
mill lattice with qqq symmetry is found to admit only 
two types of magnetic spin substructures: one of which 
can repeat as it is from one unit cell to the next, and an- 



other which undergoes a 120° spin rotation in progress- 
ing by one unit cell along the (111) direction. Thus one 
has topologically distinct ground state sectors, unit cell 
structures which are shown in Table |TT] and labeled Sq=o 
and S5=27r/3 respectively, where the ordering wavevector 

in Sg^27r/3 is Qsp — (if ; T"! T")- Curiously, as each spin 
structure features coplanar 120° rotated spins on each tri- 
angle, generic ground state spin structures are found to 
interchange between these two orderings. The only con- 
straint arising for a Sq=27r/3 spin structure is that the spin 
structure must complete three times an integral number 
of (7 = ^ unit cells within its boundaries. There are a 
macroscopic number of degenerate ground states which 
can be generated by this mixing of these two states, all of 
which exhibit a coplanar order. This tendency to form a 
coplanar spin structure is reminiscent of the nematically 
ordered state recently found for the classical Heisenberg 
model on the hyperkagome lattice.— We estimate that 
the degeneracy of the lattice grows exponentially in the 
linear lattice size, L, roughly as e^'^^^/G, as shown in 
Appendix A. 
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TABLE II: Connections and ground state candidates of the 
classical Heisenberg model on the distorted windmill (/9-Mn) 

lattice. When the lattice parameter, y = ^— j^^, then the 
nearest neighbors 1 and 2 are at equivalent distances. The fi- 
nal column shows the two spin structures minimization finds. 
Here {a, /3, 7} refer to 120° rotated spins. In moving to the 
next unit cell, the s^^o structure is unchanged, while the 
Sq=27r/3 structure replaces q — > /3 — > 7 — > a. 

Let us refer to spin configurations which order with 
ordering wavevectors, Qo and Qsp, as Qo and Qsp states 
respectively. To understand the selection of the Qo state 
over other states, we carried out MC simulations with a 
Qsp state and a Qmix state that has three planes with 
Sq={) spin configurations and three planes with Sq^2ir/3 
spin configurations for L = 6 as initial states, which al- 
lows the system to remain in these states. We found that 
the energy of the Qsp state just below the transition tem- 
perature is lower than the energy of the Qmix state, and 
the energy of the Qmix state is lower than the energy of 
Qo, as shown in Fig. [S] (a). However, the latter state is 
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FIG. 6: (Color online) For lattice size L = 6, (a) Energy and 
(b) Heat capacity vs. temperature for the AF Heisenberg 
model on the distorted windmill lattice with different initial 
spin configurations. The filled squares, circles, and triangles 
represent the Qo, Qsp, and Qmix spin configurations, respec- 
tively. The system remains in its initial configuration below 
the transition. Note that starting with any random spin con- 
figurations leads to the Qo state, (a) shows that the Qsp and 
Qmix states have the same energy as the Qo state as T — > 0, 
and have lower energy than the Qo state just below the tran- 
sition temperature. However, (b) indicates that the Qsp state 
has lower entropy below the transition than the Qmix state 
and that state in turn has lower entropy than Qo state. This 
is evidence that the selection of the Qo state is of entropic 
origin, i.e., order by disorder. 



always selected as a ground state if an initial state has 
a random configuration. This implies that the selection 
between Qo and other states is due to entropy. To con- 
firm our intuition, we computed the specific heats of the 
Qoj Qmix, and Qsp states, and found that the specific 
heat of Qo state is larger than that of Qmix or Qsp states 
below the transition temperature as shown in Fig. [S] (b). 
This confirms there are more low lying modes for the Qo 
state, which leads to the selection of the Qo state over 
the other states of the degenerate ground state manifold. 

Now let us study the cooperative paramagnetic state 
above the transition temperature. Magnetic suscepti- 
bility (not shown) has been calculated within both ap- 
proaches and shows good agreement until very close to 
the ordering temperature ~ 0.155J, well below the 
Curie-Weiss temperature Ocw — 2J, and in the zero 
temperature limit. Qualitative comparisons between the 
large- iV theory and Monte Carlo results for the structure 
factor are shown in Fig. [7] in the cooperative paramag- 
netic phase at temperatures slightly above the onset of 




FIG. 7: Contour plots of the intensity of the structure factor 
in the hhk ((a) and (b)) and hkO ((c) and (d)) planes for the 
/3-Mn lattice in the cooperative paramagnetic regime (for the 
lattice parameter y = {9 ~ ^/33)/16, L = 8 and T = 0.2J). 
Classical Monte Carlo ((a) and (c)) agrees well with large- 
TV results ((b) and (d)). The maximal intensity is shown in 
white. Axes run from (— 47r,47r) in k and h, where A;(001) 
describes the vertical axis of (a) and (b), and ft(llO) the hor- 
izontal; fc(OlO) describes the vertical axis of (c) and (d) and 
/i(lOO) the horizontal. 



magnetic order. We see excellent qualitative agreement 
between these two approaches. On the distorted windmill 
lattice the structure factor intensity is concentrated along 
the degenerate lines of the large- result. Within both 
approaches there is very little weight in the first Brillouin 
zone due to the geometric factor. In Fig. [8] a quantita- 
tive comparison along several high symmetry directions 
between these two approaches is presented. There is ex- 
tremely good agreement at all wavevectors. 



V. DISCUSSION AND SUMMARY 

A magnetic material where the magnetic properties 
can be described by the AF Heisenberg model on the tril- 
lium lattice is CelrSi. It shows^^ a Curie- Weiss suscepti- 
bility above 100 K with a magnetic moment of 2.56/.ts/Ce 
atom and a 9cw = -24 K. At low temperatures, 
shows a gradual downturn, ^^"^^ a characteristic common 
to many triangle-based frustrated magnets. X-ray scat- 
tering shows the lattice structure of this material to have 
a cubic P2i3 symmetrj*^^, with a positional parameter 
a;ce=0.6183i^ Then the Ce sites, which hold magnetic 
moments consistent with 4/^ electrons, form a trillium 
lattice of corner-shared triangles. The (non-zeroed) in- 



8 




2jl 471 671 871 1071 1271 1471 1671 
9[100] 




271 471 67: 871 1071 1271 1471 1671 

qum 




271 471 67: 871 1071 1271 1471 1671 



FIG. 8: (Color online) A quantitative comparison of the angle- 
dependent structure factor is shown along three high symme- 
try directions between the Monte Carlo (Heisenberg model) 
and large- A'^ results for classical spins in the cooperative para- 
magnetic phase on /3-Mn lattice (L = 8). The solid line is the 
large- A'^ result, while the symbols are from MC simulations. 



teraction matrix is identical to that presented in Eq. 5 of 
Ref. 13 if one takes the lattice parameter uce=xce — ^-^^ 
Curiously the energetics of the AF Heisenberg model do 
not depend on the lattice parameter u directly. As shown 
in Fig. [21 one expects that the reahzation of this model 
with different values of u will again access a cooperative 
paramagnetic regime, but with differing weights over the 
surface of the MF degenerate spheres. 

To our knowledge Ref. [2^ is the only study of the 
physics of CelrSi. In this work the magnetic susceptibil- 
ity and x-ray spectra of polycrystalline powdered samples 
have been measured, giving evidence that this material 
remains disordered to low temperatures. It would be 
very interesting to see whether or not neutron scattering 
measurement on a single crystal of CelrSi may show evi- 
dence of partial order at low temperatures giving way to 
long range order, as we predict for an antiferromagnetic 
Heisenberg model on the trillium lattice. 



The relevance of our study to /3-Mn and MnSi is 
less obvious, as both materials are metallic. However, 
the magnetic sites of each lattice features one of the 
three-dimensional corner-shared triangle lattices struc- 
tures here studied, with /3-Mn being the only known 
material to form in the distorted windmill structure. 
While the origin of the unusual non-Fermi liquid resis- 
tivity seen under pressure in MnSi^^ (Ap w AT^) is 
not at present understood, the same temperature expo- 
nent is observec^ in /3-Mn, where it is expected to result 
from antiferromagnetic spin fluctuations. We believe it 
is worthwhile to investigate the possible link between the 
magnetic fluctuations and the non-Fermi liquid behavior 
in these materials. In this vein, it is interesting to note 
that powder neutron scattering down to 1.4 K— in /3-Mn 
shows no signature of magnetic order. We are hopeful 
that this study might provide the motivation for single 
crystal neutron scattering to be carried out on (3-Mn. 

In summary, we have used large- iV theory for 0{N) 
vector spins and classical MC simulations to study the 
AF Heisenberg model on two three-dimensional corner- 
shared triangle lattices, each site of which belongs to 
three equilateral triangles. The large-A^ studies sug- 
gested that the geometrical frustration present would 
lead to a partially ordered state on both lattices. How- 
ever, through the minimization of finite size spin clus- 
ters, we found the ground state manifolds on these two 
lattices to be quite different, despite the local similarity 
between these corner sharing triangle lattice structures. 
In both cases, we found that there is a first order tran- 
sition to a magnetically ordered state using MC meth- 
ods. We further showed that the trillium lattice exhibits 
a unique ground state with a spiral ordering, while the 
distorted windmill has a macroscopic ground state degen- 
eracy. Magnetic ordering of the classical AF Heisenberg 
model on the distorted windmill lattice is therefore seen 
to arise via an order by disorder mechanism. The degen- 
eracy of this model on the trillium lattice is seen to be 
an artificial effect of the soft-constraint of the large-A^ 
theory. 

Despite the above noted differences at low tempera- 
tures between large- and MC results, in the cooperative 
paramagnetic phase above the transition temperature, 
we find a remarkable resemblance between the respec- 
tive spin-spin correlations. This leads us to ask whether 
the salient features of the large-iV theory, the angular 
and directional dependences of the spin-spin correlations 
found in the partially ordered state obtained by large- iV 
theory are present at finite temperatures above the tran- 
sition temperature. As true partial order exhibits long 
range order along particular directions only as T — s- 0, 
it is not possible to have partial order at any finite tem- 
peratures, since the spin-spin correlation decay exponen- 
tially at finite temperatures. This being said, the quali- 
tative directional dependence characteristic of a partially 
ordered state survives above the transition temperature 
(note that the ground state is smoothly connected to the 
cooperative paramagnetic phase in large-A'^ theory) al- 
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FIG. 9: Contour plots of the intensity of the large-A'^ theory structure factor in the hhk plane for different values of u. The 
vertical axis is fc(OOl), and the horizontal axis is ft(llO), where k and h both run from (— 47r,47r). Varying the parameter u, 
does not change the energy within either the large- theory (here N — 3, L — 24 and T = 1/4. 5 J), nor in classical Monte 
Carlo simulations. It does, however, change the geometric contribution to the structure factor quite dramatically. While near 
u = 0.125 (center) one appears to see few vestiges of the degenerate spheres of the treatment, for both large (right) and small 
(left) values of u, one sees that these are quite well captured. Curiously, most real systems seem to lie closer to the middle 
graph. 



lowing us to conclude that a "disguised" partial order has 
been recovered in the cooperative paramagnetic phase. 
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APPENDIX A: DEGENERACY OF THE 
DISTORTED WINDMILL LATTICE GROUND 
STATE 

To estimate the degeneracy of the ground state, we re- 
alize that fixing {a,(3,j} in Table H] leaves us with six 
tilings of the spins in the unit cell, which we might la- 
bel {Sg=0' 'Sq^OJ ■Sg=OJ ■Sg=27r/3' '^?=2Tr/3' '^?=2Tr/3}- ^^^^ * 

is taken to be the spin structure as presented in Table HI 
the same taking a — > /3 — > 7 — s- a, and the same 
taking a j f3 a. In progressing from one unit cell 
to the next along the (111) direction, we can follow s^^q 
with cither sJ^q or sJ^2,r/3' ^q=2n/3 '^^^ followed by 
either s^^27r/3 ^q=o- Replacing 1 ^ 2 ^ 3 ^ 1 gives 
the general rules for allowed spin structures (provided 
the number of Sg=27r/3 structures is a multiple of three). 



With these rules, it becomes a problem in combinatorics 
to determine the ground state degeneracy. 

For simplicity, '^^ let us consider a finite size lattice with 
LxLxL unit cells, with L divisible by 3. Clearly we have 
one state with only Sq^27r/3 ^P^"^ structures. By remov- 
ing three s^^2Tr/3 ^P^^ structures, we can insert three s;^^q 
spin structures in their place, or more generally, remov- 
ing 3m planes of unit cells which rotate as they progress 
along the (111) direction, allows the introduction of 3m 
planes of spin structures which need not rotate from one 
plane to the next. Upon such a removal there should be 
31 — 3m available unit cell planes for 3m equivalent spin 
structures. This implies a degeneracy of approximately^ 



(3Z-1)! 



(3m)!(3Z-3m- 1)! 



(Al) 



where m can be {0, 1, .., I}. Summing over m then gives 
the total number of states: 



Nr 



E 



(3/-1)! 



(3m)!(3;-3m-l)! 



3/-1 

E 



(3/-1)! 



23/- 



3 ^ m!(3/ - m - 1)! 

m=0 ^ ' 
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(A2) 



(A3) 



so the number of ground states grows exponentially m 

. . „ln2I, 0.691, 

the hnear lattice size, roughly as A^ « « ^-j^ — . 
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In the graphs here we have used it = 0.138 which 
is the value relevant for the Mn sites on the 
MnSi lattice. This corresponds to the distance 

^(i)2-H(i-2u)2-h(2u)2. Taking uce = 0.1183 

describes the nearest neighbour distances between 



Ce spins since 



y^(i)2 + (I - 2k)2 + (1 - 2«)2 = 



y^(i)2 + (|-2(x - i)2 + {2{x - i))2. 

In the thermodynamic limit corrections resulting from this 
assumption will be small. 

Some of these states may be equivalent, but we expect 



this to contribute at most a factor of ■ 



31 — 3m 

merely be a logarithmic correction to our result 



which should 



